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Abstract 

Using the SU(2) lattice QCD, we formulate the dual Wilson loop and study 
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theory including the electric current j„ and the monopole current k^. After 
the abelian projection in the MA gauge, the system can be separated into 
the photon part and the monopole part corresponding to the separation of j M 
and k^, respectively. We study here the monopole part (the monopole-current 
system), which is responsible to the electric confinement. Owing to the ab- 
sence of electric currents, the monopole part is naturally described using the 
dual gluon field without the Dirac-string singularity. Defining the dual 
Wilson loop from the dual gluon B^, we find the perimeter law of the dual 
Wilson loop in the lattice QCD simulation. In the monopole part in the MA 
gauge, the inter-monopole potential is found to be flat, and can be fitted as 
the Yukawa potential in the infrared region after the subtraction of the arti- 
ficial finite-size effect on the dual Wilson loop. Prom more detailed analysis 
of the inter-monopole potential considering the monopole size, we estimate 
the effective dual-gluon mass ms — 0.5GeV and the effective monopole size 
R ~ 0.2fm. The effective mass of the dual gluon field at the long distance can 



be regarded as an evidence of "infrared monopole condensation" . 
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I. INTRODUCTION 



Quantum chromo-dynamics (QCD) is the fundamental theory of the strong interaction 
and is an SU(iVc) nonabelian gauge theory described by the quark field q and the 
gluon field as 

£qcd = —tviG^Gn + m - rn q )q, (1) 

where is the SU(iV c ) field strength = j^[D^,D u ] with the covariant-derivative 
operator = d^+ieA^. Due to the asymptotic freedom, which is one of the most important 
features in QCD the gauge-coupling constant of QCD becomes small in the ultraviolet region 
Accordingly, the perturbative QCD can describe the high-energy phenomena like the 
Bjorken scaling and the hadron jet properties 

On the other hand, in the low-energy region, the QCD-coupling constant becomes large, 
and there arise the nonperturbative-QCD (NP-QCD) phenomena such as color confine- 
ment and dynamical chiral-symmetry breaking corresponding to the strong-coupling nature. 
These NP-QCD phenomena are extremely difficult to understand in the analytical manner 
from QCD, and have been studied by using the effective models or the lattice QCD sim- 
ulation [^]. Here, the lattice QCD Monte Carlo simulation is the numerical calculation of 
the QCD partition functional, and it is one of the most reliable methods directly based on 
QCD. In fact, the lattice QCD simulations well reproduce nonperturbative quantities such 
as the quark static potential, the chiral condensate (qq) and low- lying hadron masses @. 

Recently, the lattice QCD simulation has shed light on the confinement mechanism in 
terms of the dual-superconductor picture, which was proposed by Nambu,'t Hooft and Man- 
delstam in the middle of 1970's In this scenario, quark confinement can be understood 

with the dual version of the superconductivity. In the ordinary superconductor, the Meiss- 
ner effect occurs by condensation of the Cooper-pair with the electric charge. Consider 
the existence of the magnetic charges with the opposite sign immersed in the superconduc- 
tor, then the magnetic flux is squeezed like a tube between the magnetic charges, and the 
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magnetic potential between them becomes linear as the result of the Meissner effect ITT . 
In the dual-superconductor scenario, the QCD vacuum is assumed as the dual version of 
the superconductor, and the dual Meissner effect brings the one- dimensional flux squeezing 
fT2|| between the quark and the anti-quark, which leads to the linear confinement potential 

H0 



The dual Higgs mechanism, however, requires "color-magnetic monopole condensation" 
as the dual version of electric condensation in the superconductor, although QCD dose 
not include the color-magnetic monopole as the elementary degrees of freedom. On the 
appearance of magnetic monopoles from QCD, 't Hooft showed that QCD is reduced to an 
abelian gauge theory with magnetic monopoles by taking the abelian gauge, which fixes the 
partial gauge symmetry SU(iVc) /U(l) Arc_1 through the diagonalization of a gauge-dependent 
variable |14[- Here, the monopole appears as the topological object corresponding to the 



nontrivial homotopy group n2(SU(N e )/U(i) Na - 1 )=Z** e - 1 . 

As for the irrelevance of off-diagonal gluons, recent lattice QCD studies show the abelian 
dominance \T5\ for the NP-QCD phenomena in the maximally abelian (MA) gauge |TB|,|T7 



For instance, confinement P^JT9[| and dynamical chiral-symmetry breaking pO| , pT| are almost 
described only by the diagonal gluon component, in the MA gauge. Then, taking the MA 
gauge and removing off-diagonal gluons, the abelian-projected QCD (AP-QCD) is obtained 
as the abelian gauge theory keeping the NP-QCD features. AP-QCD includes not only 
the electric current but also the magnetic current k^, and can be decomposed into the 
monopole part and the photon part corresponding to the separation of and j^, respectively 
p3fl. The lattice QCD studies show that only the monopole part is responsible to NP- 



QCD phenomena pm|23| -r2q1 especially to the electric confinement p3|-f25[] in the MA gauge. 
This is called as the monopole dominance. 

In this paper, we concentrate the monopole part (the monopole-current system) in the 
MA gauge in QCD, and study the dual Higgs mechanism in the QCD vacuum based on the 
dual gauge formalism pTySSfl . To this end, we perform the SU(2) lattice QCD simulation 
in the MA gauge, and extract the monopole current k^ as the relevant degrees of freedom 



for electric confinement. Then, we calculate inter- mo nopole potential in the monopole part, 
(the monopole current system) in QCD to examine monopole condensation, and evaluate 
the effective mass of the dual gluon field p7|j28[1 . 



II. SEPARATION OF AP-QCD INTO THE MONOPOLE PART AND THE 
PHOTON PART IN THE MA GAUGE 

A. MA Gauge Fixing and AP-QCD 

Recent studies with the lattice QCD Monte Carlo simulation have revealed the abelian 
dominance and the monopole dominance in the maximally abelian (MA) gauge for the non- 
perturbative QCD (NP-QCD) phenomena such as confinement, dynamical chiral-symmetry 
breaking and instantons [18-26]. 

In the continuum Euclidean QCD with N c = 2, the MA gauge fixing is defined by the 
minimizing the total amount of off-diagonal gluons, 

R oS = J d*x[{A»{x)} 2 + {A»{x)} 2 ] = 2 J d 4 x[Al(x)A^x)] (2) 

with 

M x ) =y k A »T> (3) 

a=l Z 

A$(x) = ±{Al{x)±iAl(x)}, (4) 

by the SU(2)-gauge transformation. In the MA gauge, off-diagonal gluon components, A 1 
and A 2 , become as small as possible, and the gluon field A^ mostly approaches to the abelian 
gluon field, A^ hel = A^y, In the MA gauge, the SU(2) local symmetry is reduced into the 
U(l)3 local symmetry with the global Weyl symmetry. Under the residual C/(l)3-gauge 
transformation with u = exp^—icj) 1 ^-) £ U(lj3, the gluon components are transformed as 

Al(x)^Al(x) + d^(x), (5) 
A±{x) - eW*)A$(x). (6) 
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Then, in the MA gauge, the diagonal gluon A 3 ^ behaves as the abelian gauge field, while off- 
diagonal gluons behave as charged matter fields in terms of the residual gauge symmetry. 

As a remarkable feature of the MA gauge, the abelian dominance holds for the NP- 
QCD phenomena such as quark confinement and chiral-symmetry breaking [18-21]. Here, 
we call abelian dominance for an operator O, when the expectation value (OL4 M ]) is almost 
equal to the expectation value {0[A Ahel ]) M \, where off-diagonal gluons are dropped off in 
the MA gauge. For instance, the abelian string tension oam = (a(A Ah el )) MA in the MA 
gauge is almost equal to <7su(2) = (°"(^ U( ' 2 ' ) )) as OAbci — 0.92<t S u(2) for (3 ~ 2.5 in the lattice 
QCD [24,25]. Thus, NP-QCD phenomena are almost reproduced only by the abelian gluon 
A^ hcl , and off-diagonal gluon components A^ do not contribute to NP-QCD in the MA 
gauge. Hence, as long as the infrared physics is concerned, QCD in the MA gauge can be 
approximated by the abelian projected QCD (AP-QCD), where the SU(2) gluon field A^ 
is replaced by the abelian gluon field A Abel . In other word, AP-QCD is the abelian gauge 
theory keeping essence of NP-QCD, and is extracted from QCD in the MA gauge. Hereafter, 
we pay attention to the AP-QCD described by A Ahel in the MA gauge. 

The abelian- projected QCD (AP-QCD) includes not only j M but also k^. In this subsec- 
tion, we investigate the general argument on the extended electro-magnetic system including 
both the electric current j M and the magnetic current k^. In the extended Maxwell equations 
with jn and k^ the field strength F^ el satisfies as 

d,F^ = j v (7) 

% *i^ bcl = ~d, V = K (8) 

with *F^, bel = ^nvpoF^o^ '• In the presence of both and k^, the field strength F^ el cannot 
be described by the simple two-form d^A^ 1 — d u A Ahel with the regular one-form A Ahel [29]. 
In general, the field strength F^ el consists of two parts, 

F^ cl =(dAA Ahd ), u + V, (9) 

where the former part denotes the ordinary two-form and the latter part = ^e^a/sCap 
denotes the Dirac-string singularity [29]. Here, £^ can be written as 
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^u(x) = — -(n A k),,,, 



1 

n ■ d 

1 

n ■ d 



d4 y( x \z—5\y)( n v k »(y) - ^k^y)) 



= J d 4 y9(x n - 2/ n )(5 3 (x x - y L )(nMy) ~ n*k^y)), (10) 

where x n = x^n^ and xj^ = x M — (x ■ riju^. Here, is arbitrary four- dimensional unit 
vector corresponding to the direction of the Dirac string. Thus, in the ordinary description 
[29] the system includes the singularity as the Dirac string which makes the analysis 
complicated. 



B. Separation of AP-QCD into Photon Part and Monopole Part 

To clarify the roles of and k^ to the nonperturbative quantities of QCD, we consider 
the decomposition of AP-QCD into the photon part and the monopole part, corresponding 
to the separation of j M and k^. We call this separation into the photon and monopole parts 
as the "photon projection" and the "monopole projection", respectively. [26,30]. 

The field strength F^ el in AP-QCD is separated into in the photon part and F^° 
in the monopole part 

F^ cl = + F$>. (11) 
As the physical requirement, Fjj 1 and F^f satisfies the Maxwell equations as 

d lt F$ = j v , d,*F^ h = 0, (12) 
«V^° = 0, dp*F*> = k u , (13) 

respectively. From the requirement (jl2|), there exists the regular vector A^ h is defined so as 
to satisfying 

= d,Al h - <9,< h (14) 
in the photon part. This relation is the same as F^ = d ll A u — d v A in the ordinary QED. 
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As for the monopole part, cannot be expressed by the simple two- form, but is 
expressed as 

d»A™° - d v Af° = F™° + V> (15) 

where the second term provides the breaking of the Bianchi identity. From Eqs.(§), (|lT|), 
(0) and (|15D , one finds the relation satisfy 



d A A Ahcl = d A (A Fh + A Mo ), (16) 

and then the difference between A Ahel and A^+A^ 10 is at most the total differential term 
as d^x- Hence, we can set 

Af eI = A™ + A™° (17) 

by taking a suitable gauge without loss of generality. 

In the practical calculation, the photon part variable A Fh and the monopole part variable 
A^}° can be obtained as 



Al h = u^j v = a-%F^ e \ (18) 



A Mo = n -lQ ^ (1Q) 



where the inverse d'Alembertian □ 1 is the nonlocal operator [31], 



<*I D_1 W = (20) 



which satisfies D x (x\0~ 1 \y) = (x\y) = 5 4 (x — y). 

Let us check that A Ph and A^° defined in Eqs.(|TB]) and (|T^P, satisfies the physical re- 
quirement in Eqs.(14) and (15) First, we consider the photon part with A Fh . Starting from 



Eq.(18), one easily finds d^A^ 1 = 0, and hence the first Maxwell equation in Eq. ([T2|) in the 



photon part can be derived as 



d,F™ = d,(d,A p u h - d u Al h ) = DA™ - d u (d,Al h ) = UA? 

= UU- 1 j v =j v . (21) 



The second Maxwell equation in Eq. flT2|) is automatically derived as the Bianchi identity for 
the two form of A^ h . Thus, the photon part with A^ h defined by Eq.(Jf^) does not include 
the magnetic current but only includes the electric current. 

Second, we consider the monopole part with A^°. From Eqs.(|9|), (|T8| ) and (|T9|), the sum 
of A^, h and can be written as 

A Ph + A Mo = D-l^Abel + = n -l da{daA AM _ ^Abel) 

= A^ hEl -d,(d a At hel ) } (22) 
and hence the two-form of A^}° satisfies the physical requirement (|15|) as 
{d^° - d v Af°) = d,At bcl - d v Af>« 1 - (dX h ~ d»Al h ) 

= ^+*U- (23) 



Therefore, in the monopole part, the Maxwell equations is derived as 

d»F™° = d,F^ cl - d»F™ = 0, (24) 
d, *F™° = d, *i^ bel - d, *F™ = k v , (25) 

starting from Eq. ([191). 

Thus, A^ 1 and A^° defined as Eqs.(|l8D, and ( |i~9D satisfy the physical requirement 
Eqs.(|T4]) and (p!5[). Hence, the monopole part carries the same amount of the magnetic 
current as that is the original abelian sector, whereas it dose not carries the electric current. 
The situation is just the opposite in the photon part. In the actual lattice QCD simula- 
tion, the monopole current and the electric current are slightly modified through the 
monopole and the photon projections, respectively, due to the numerical error on the lattice. 
However, these differences are negligibly small in the actual lattice QCD simulation. In fact, 
k^f° ~ and ~ hold in the monopole part, and ~ j M and k^ h ~ hold in the 
photon part within 1% error. Here, we have kept the labels as "Mo" and "P/t" for the 
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electric current and the monopole current, and we have used (k^°,jjf°) and (k^ h ,j^ h ) for 
these currents in the monopole part and the photon part, respectively. 

As a remarkable fact, lattice QCD simulations show that nonperturbative quantities such 
as the string tension, the chiral condensate and instantons are almost reproduced only by 
the monopole part in the MA gauge, which is called as monopole dominance [23-26]. On the 
other hand, the photon part dose not contribute these nonperturbative quantities in QCD. 

Since we are interested in the NP-QCD phenomena, it is convenient and transparent 
to extract the relevant degrees of freedom for NP-QCD by removing irrelevant degrees of 
freedom like the off-diagonal gluons and the electric current j M . Therefore, we concentrate 
ourselves to the monopole part, which keeps the essence of NP-QCD as confinement. 

III. DUAL GAUGE FORMALISM 
-DUAL GLUON FIELD AND DUAL WILSON LOOP- 

In this section, we study the monopole part of the QCD vacuum using the dual gauge 
formalism [27,28]. In the MA gauge, the monopole part carries essence of the nonpertur- 
bative QCD as the electric confinement. According to the absence of the electric current 
(jfi — 0), the Maxwell equation in the monopole part becomes 

d»F$> = (26) 
= k u , (27) 

where denotes the field strength in the monopole part. This system resembles the dual 
version of QED with ^ and k^ = 0, and hence it is useful to introduce the dual gluon 
field in the monopole part for the study of the dual Higgs mechanism in QCD [27,28]. 
The dual gluon field is defined so as to satisfy the relation 

d,B v -d v B^*F^\ (28) 

which is the dual version of the ordinary relation F^ = - d v Af° in QED. The 

interchange between and corresponds to the electro-magnetic duality transformation, 
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F^v <-> *F ilv or H <-> E. Owing to the absence of j M , the dual gauge field 5 M can be 
introduced without the Dirac-string singularity. In the other words, the absence of j u is 
automatically derived as the dual Bianchi identity, 

Ju = d,F, u = d,\dAB), u = 0. (29) 

Let us consider the derivation of the dual gauge field _B M from the monopole current 
k^. Taking the dual Landau gauge d^B^ = 0, the Maxwell equation *F^° = d 2 B v — 
dvid^By) = k v is simply reduced as UB^ = k^. Therefore, the dual gluon field B^ is 
obtained by using the inverse d'Alembertian as 

B v (x) = U- l k u (30) 

or equivalently 

B v (x) = J d 4 y(x\n- l \y)k„(y) = ~J d 4 yj-±-^h(y). (31) 

Thus, the monopole part is described by the monopole current /c M and the dual gluon £? M in 
the regular manner based on the dual gauge formalism. 

In the dual superconductor picture in QCD, k^ and B^ correspond to the Cooper-pair 
and the photon in the superconductor, respectively. The Cooper-pair and the photon are 
essential degrees of freedom which bring the superconductivity. In the superconductor, the 
photon field gets the effective mass as the result of Cooper-pair condensation, and this 
leads to the Meissner effect. Accordingly, the potential between the static electric charges 
becomes the Yukawa potential Vy(r) oc in the ideal superconductor obeying the London 
equation. Similarly, the dual gluon B^ is expected to be massive in the the monopole- 
condensed system, and the mass acquirement of B^ leads to the dual Meissner effect. In 
other words, the acquirement of dual gluon mass m B reflects monopole condensation, and 
brings electric confinement. Hence, we can investigate the dual Higgs mechanism in QCD by 
evaluating the dual gluon mass tub, which is estimated from the inter-monopole potential. 
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To estimate the interaction between the monopoles, we propose the dual Wilson loop 
W D [27,28]. The dual Wilson loop W D is defined by the line- integral of the dual gluon field 
B^ = B^- along a closed loop C, 

W D (C) = hi exp (te M j c B^dx^) = Re[exp (i^ £ Bjtte^)], (32) 
which is the o?ua/ version of the abelian Wilson loop 

W Ahel (C) = ^trexp(ie£ cl rf^) = i?e[exp (i| £ Ajd^)]. (33) 

Using the Stokes theorem, the dual Wilson loop Wd(C) is rewritten as 

W D (C) = itrexp (ie M £ *F™°dS, v ), (34) 

with the dual gauge field strength Fj^,°. The dual Wilson loop Wd(R x T) describes the 
interaction between the monopole-pair with the test magnetic charges ^f- and — ^f-. Here, 
these test magnetic charges are pair-created at t — and are pair-annihilated at t — T 
keeping the spatial distance R for < t < T. As T goes to infinity, the dual Wilson loop 
(W D (R x T)) means the interaction between the static monopole and anti-monopole with 
the separation of the distance R. The inter-monopole potential is obtained from the dual 
Wilson loop as 

V M (R) = ~ Jim i \n(W D (R x T)> (35) 

in a similar manner to the extraction of the inter-quark potential from the Wilson loop 
[4-7]. To summarize here, for the investigation of the dual Higgs mechanism in QCD, we 
have introduced the dual gluon field B^ and the dual Wilson loop in the monopole part of 
the AP-QCD in the MA gauge. In the next section, we consider the practical procedure on 
the calculation of the dual Wilson loop and the inter-monopole potential in the lattice QCD 
formalism. 
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IV. DUAL GAUGE FORMALISM ON THE LATTICE 



We study the dual Wilson loop and the inter-monopole potential in the MA gauge using 
the SU(2) lattice QCD Monte Carlo simulation. The lattice QCD simulation is the direct 
calculation of the QCD partition functional using the Monte Carlo method. The physical 
expectation value of an observable is numerically obtained by averaging its value over 

all gauge configurations with the weight factor exp(— SfA^]), where S[Afj] denotes the lattice 
QCD action in the Euclidean metric J7[. 

In the lattice gauge formalism with the lattice spacing a, the SU(2) link variable is 
defined by U^(s) = exp(iaeA^(s)) = exp(meyl^(s)r a /2) G SU(2), where e and r a /2 denote 
the QCD gauge coupling and the generator of the SU(2) group, respectively. The standard 
lattice QCD action in the gauge sector is defined by 

^^Eli-^M + Wl' P=^T ( 36 ) 

s,n,v /iV c e 

using the plaquette variable U^ u (s) = U^sjU^s + //)L^(s + u)UUs). In the continuum 
limit a — > 0, the plaquette U^ v (s) becomes exp[ia 2 G^ u (s)}, and hence S L coincides with the 
continuum QCD action 

S = J d A xUi{G^G, v ). (37) 

Thus, the QCD system is described by the link variable U^s) G SU(iV c ) instead of the gauge 
field A^(x) G su(N c ) in the lattice formalism. 

Here, we consider the extraction of the abelian-projected QCD (AP-QCD) from the 
lattice QCD. In the SU(2) lattice formalism, the MA gauge fixing is achieved by maximizing 



R = hrj:[r 3 U,(s)r 3 Ul(s)] ="£[1-2 ({£/»} 2 + {t/»} 2 )] (38) 
by the SU(2) gauge transformation, 

U^Uf A = V(s)U,{s)V\s + fi) } (39) 
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where V(s) and V(s + fi) are the gauge functions located at the starting and end points of 
the link variable U^(s). In this gauge, the absolute value of off-diagonal components U^(s) 
and £7*(s) are forced to be small as possible using the gauge degrees of freedom. 

In accordance with the Cartan decomposition, the SU(2) link variable U^s) is factorized 



as 



U, 



MA, 



M„ 



exp^r 1 + 9lr% u,(s) ee expMr 3 )] 



(40) 
(41) 



where u M £ U(l) 3 and M M £ SU(2)/U(1) 3 correspond to the diagonal part and the off- 
diagonal part of the gluon field, respectively. In the continuum limit, the angle variable d a 



goes to the gluon field A a as Q a — > ^eaA a . The off-diagonal factor M M (s) is rewritten as 



with 



^ cos 9,, — sin^,, e * Xm ^ 
sin6 l „e* XM cos#„ 



(42) 



r 7T , 



X = tan . 



(43) 



Under the abelian gauge transformation with v(s) £ U(l)3, M M (s) and u^(s) are transformed 

as 



M fl (s)^M;{s) = v{s)M ll ( S )v\s), 



u.As: 



K(s) = v(s)u fl (s)v*(s + p), 



(44) 
(45) 



to keep the form of Eq.(gl]) for £ SU(2)/U(1) 3 and £ U(l) 3 . Then, M^s) behaves 
as the charged matter field and the abelian link- variable 



I iff 3 n \ 

e c 



V 



e" 



(46) 



behaves as a abelian gauge field with respect to the residual abelian gauge symmetry. 
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In the lattice QCD, the abelian dominance is expressed as (O(L^)) ~ (0(m m ))ma for the 
infrared quantities such as a and (qq) in the MA gauge. The abelian projection is performed 
by the replacement of U^(s) — > u^s), and then the abelian-projected QCD (AP-QCD) on 
the lattice is described only with the abelian link- variable which is the diagonal part 

of the SU(2) link- variable Ujf A (s) in the MA gauge. 

Next, we consider the separation of AP-QCD into the monopole part and the photon 
part in the lattice formalism. Using the diagonal gluon component 0^(s), the abelian field 
strength 9^ el is defined by 

0^ 1 = mod 27r (^-<9,0j) (47) 
= dX(s) ~ d v 9l(s) + 27T7v(s), (48) 

where the former part denotes the ordinary two-form and n^s) G Z corresponds to the 
Dirac string on the lattice 123]. In the lattice formalism, the photon part 0„ h {s) and the 



monopole part 0^ Io (s) are obtained from 9^ cl (s) and 2irn^ u (s), respectively 

<%\s) = -{□- 1 ^C Cl }(*) (49) 
6™ {s) = -2n{a~ l d„n^}(s), (50) 

using the inverse d'Alembertian D _1 on the lattice [22-28]: The diagonal gluon component 
6^(s) is found to be decomposed as 

ol(s) = e^(s) + e^(s) (si) 

in the Landau gauge, d^Js) = 0. 

The field strengths, 8^ in the photon part and in the monopole part are given as 

0™ = mod 2n (dX h ~ dj™) (52) 
Op = mod 27r (9 M C° - 9*0?°) (53) 

on the lattice. In the continuum limit a — > 0, these field strengths becomes as 
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)Mo 

jj,u 



)Ph 




(54) 



(55) 



Using the dual field strength *9 



= je^ upa 9pa°, the dual Wilson loop is expressed as 



W D (C) = ^trexp (i [ 



\ Js 




) 



(56) 



using the Stokes theorem on the lattice. 

To summarize, we show the procedure on the derivation of the dual Wilson loop from 
lattice QCD as follows: 

1. We generate the SU(2) gauge configurations {?7 M (s)}i using the Monte Carlo method 
for the lattice QCD. 

2. We carry out the gauge transformation, U^(s) — > Ujf A (s), so as to satisfy the MA 
gauge fixing condition with the minimization of R. 

3. The SU(2) link-variable U^s) is factorized as U^s) = M fl (s)u IJi (s), and the abelian 
projection is performed by the replacement of U^(s) G SU(2) by the abelian link- 
variable Un(s) = exp{z6^(s)r 3 } G U(l)3. 

4. The two-form of the diagonal gluon component is decomposed as dffi(s) — d v 9^(s) = 
9f^ el (s) +2'im lu ,(s) with the abelian field strength 9^ el G (— ir, n] and the Dirac string 
lim^y G 2nZ. 

5. The abelian gauge field 9^ is decomposed as 6 3 = 9^° + 9^ h with the photon part 9^ h 
and the monopole part 9^f°, which are obtained from 9^ el and n pv by way of Eqs.(f49|) 
and ( p0| ) using the inverse d'Alembert operator in the Landau gauge. 

6. Using the field strength 9^° in the monopole part, the dual Wilson loop (W / d(C))ma 
and the inter- monopole potential Vj^(r) are calculated with Eqs.flSED and Then, 
the effective dual-gluon mass m B is estimated from the inter-monopole potential Vjy(r). 
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Thus, we investigate the dual Higgs mechanism in QCD by extracting the dual Wilson loop 
(Wd(C))ma and the inter-monopole potential Vm{t) from the gauge configurations obtained 
in the lattice QCD simulations. 

V. LATTICE QCD RESULTS FOR INTER-MONOPOLE POTENTIAL AND 

DUAL GLUON MASS 

In this section, we show the numerical result of the lattice QCD simulation. For the study 
of the dual Higgs mechanism in QCD, we calculate the dual Wilson loop W D (R,T) and the 
inter-monopole potential Vm{t) in the monopole part (the monopole- current system) in QCD 
in the MA gauge using the SU(2) lattice with 20 4 and f3 = 2.2 ~ 2.3. All measurements are 
performed at every 100 sweeps after a thermalization of 5000 sweeps using the heat-bath 
algorithm. The physical unit or the lattice spacing a is determined so as to reproduce the 
string tension a = 1 GeV/fm for each (3, e.g. a = 0.199 fm for f3 = 2.3. We prepare 
100 samples of gauge configurations. These simulations have been performed using the 
super-computer SX-4 at Osaka University. 

The dual Higgs mechanism is characterized by the effective-mass acquirement of the dual 
gluon B^, which is brought by monopole condensation. To examine the effective mass of B^, 
we calculate the inter-monopole potential V M (r) from the dual Wilson loop (W D (R,T)) MA 
obtained in the lattice QCD. As shown in Fig.l, the dual Wilson loop {Wd(R,T))ma seems 
to obey the perimeter law rather than the area law for large loops with I, J > 3. Since the 
dual Wilson loop (Wd(R,T)) M a satisfies the perimeter law as 

\n(W D (RxT)) MA ~-2(R + T)-a (57) 

for large R and T, the inter-monopole potential becomes constant 2a in the infinite limit of 
T, 

2a 

V M {R) -> lim — R + 2a = 2a. (58) 
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In the actual lattice QCD calculation, however, we have to take a finite length of T, and 
hence the linear part (2a/T)R slightly remains as a lattice artifact [28]. Therefore, we 
have to subtract this lattice artifact (2a/T)R for evaluating of the inter-monopole potential 
Vm(R) from (Wd(R,T)) in the lattice QCD simulation. Here, a can be estimated from the 
slope of the dual Wilson loop \n(W D (R x T)) MA for large R and T for each lattices. We 
obtain a = 0.141 ± 0.025GeV from the dual Wilson loop in Fig.l(b) with (3 = 2.3. 

After the subtraction of the lattice artifact (2a/T)R, we consider the slope of the inter- 
monopole potential Vm(t) in the lattice QCD. As shown in Fig. 2, the inter-monopole poten- 
tial Vm{t) is short-ranged and flat in comparison with the linear inter-quark potential with 
string tension a = lGeV/fm. Now, we try to apply the Yukawa potential Vy(r), 

to the inter-monopole potential Vm{t)- As shown in Fig. 3, the inter-monopole potential can 
be fitted by the Yukawa potential Vy(r) in the long distance region, and we evaluate the 
dual gluon mass asm B ~ 0.5GeV. 

Finally, we consider the possibility of the monopole size effect, because the QCD 
monopole is expected to be a soliton like object composed of gluons. In fact, from the 
recent lattice QCD study, the QCD monopole includes large off-diagonal gluon components 
near its center even in the MA gauge [28,30,32], and the off-diagonal gluon richness would 
provide the "effective size" of the QCD monopole similar to the 't Hooft-Polyakov monopole 
[1-4] We introduce the effective size Ru of the QCD-monopole, and assume the Gaussian- 
type distribution of the magnetic charge around its center, 



1 ,-lx 



2 



P (x;Rm)= , r- x 3 exp(— — )■ (60) 

Since the monopole part is an abelian system, simple superposition on is applicable like 

the Maxwell equation. Therefore, the inter-monopole potential with the effective monopole 

size Rm is expected to be 

t t- / D v (e/2) 2 f 3 f 3 exp(-m B |x-xi +x 2 |) 

V(x;R M ) = -. / « ^i / <rx2p(xi; i? M )p(x 2 ; R M ) ; : ; , (61) 

Air J J |x — xi + x 2 | 
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or equivalently 



V(r;R M ) 



exp[— m B J {J (r — r 2 cos 6* 2 ) 2 + (r 2 sin# 2 ) 2 — T\ cos 9\} 2 + (ri sin^i) 2 ] 




/*00 Z 1 OO P7T /*7T 

/ dn / dr 2 e- (r ? +r 2 )/fi M / ^ / d0 2 sin 0i sin 2 



(62) 



{v/(r — r 2 cos# 2 ) 2 + (r 2 sin0 2 ) 2 — ricos^i} 2 + (rising) 2 



where r = |x — y | is the distance between the two monopole centers. We apply the Yukawa- 
type potential V(r;R M ) to the inter-monopole potential Vm{t) in Fig.3. The potential 
V{j-;Rm) with the effective monopole size Rm = 0.21fm seems to fit the lattice data of 
Vm{t) in the whole region of the distance r. 

Thus, we estimate the dual gluon mass m B ^ 0.5GeV and the effective monopole size 
R M ~ 0.2fm by evaluating the inter-monopole potential Vm{t) from the dual Wilson loop 
Wd(R, T) in the monopole part in the MA gauge. In the long-distance region, we find 
the effective-mass acquirement of the dual gluon B^, which is essential for the dual Higgs 
mechanism in the dual superconductor scenario. This result suggests "infrared monopole 
condensation" or monopole condensation in the long-scale description of the QCD vacuum. 
The monopole size R M would provide a critical scale [30-32] for the nonperturbative QCD in 
terms of the dual Higgs theory, because the QCD-monopole structure such as off-diagonal 
gluons [28,30] should be considered at the shorter scale than Rm, similar to the structure of 
the 't Hooft-Polyakov monopole. 



To examine the dual superconductor picture for the quark confinement mechanism in 
the QCD vacuum, we have studied the dual Higgs mechanism in terms of the effective-mass 
acquirement of the dual gluon field using the lattice QCD Monte Carlo simulation. In 
the MA gauge, QCD is reduced to the abelian gauge theory with the color-electric current j M 
and the color- magnetic monopole current k^. The abelian-projected QCD, the diagonal part 
of QCD, can be separated into the photon part and the monopole part corresponding to the 
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separation of and k^, respectively. Reflecting the abelian dominance and the monopole 
dominance, the monopole part carries essence of NP-QCD and then is of interest in the MA 
gauge, so that we have concentrated ourselves to the monopole part (the monopole- current 
system) in QCD. 

In order to investigate the dual Higgs mechanism in QCD, we have introduced the dual 
gluon field and have studied its features in the monopole part in the MA gauge in the 
lattice QCD. Owing to the absence of the electric current, the monopole part resembles the 
dual version of QED, and hence this part is naturally described by the dual gluon field 
without meeting the difficulty on the Dirac-string singularity. In the dual gauge formalism, 
the dual Higgs mechanism is characterized by the acquirement of the effective mass of 
the dual gluon field B^. Then, to evaluate the dual gluon mass, we have calculated the dual 
Wilson loop {Wd(R x T)) M a, and have studied the inter-monopole potential Vm[t) in the 
monopole part in the MA gauge based on the dual gauge formalism by using the lattice 
QCD simulation. 

In the lattice QCD, we have found that the dual Wilson loop obeys the perimeter law for 
large loops. Considering the finite-size effect of the dual Wilson loop, we have investigated 
the inter-monopole potential Vm(V), an d have found that Vm{t) is short-ranged and flat 
in comparison with the linear inter-quark potential. Then, we have compared the inter- 
monopole potential Vm(V) with the Yukawa potential and have estimated the dual gluon 
mass as m# ~ 0.5GeV, which is consistent with the phenomenological parameter fitting 
in the dual Ginzburg-Landau theory fL3| . The generation of the dual gluon mass mg in 



the infrared region suggests the realization of the dual Higgs mechanism and monopole 
condensation in the long-scale description of the QCD vacuum. In this way, we have shown 
the evidence of "infrared monopole condensation" in the lattice QCD in the MA gauge. 

To explain the short-range deviation between the inter-monopole potential Vm{t) and 
the Yukawa potential, we have considered the effective size Ru of the monopole, since the 
monopole would be a soliton-like object composed of gluons [28,30-32]. The lattice data of 
the inter-monopole potential can be well fitted with the Yukawa-type potential V(r; Ru) 
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with the effective size Rm — 0.2 fm of the monopole. This monopole size Rm — 0.2 fm may 
provide the critical scale for the dual Higgs theory in QCD, because the monopole structure 
relating to off-diagonal gluons become visible [28,30-32] and the QCD system cannot be 
described only with the abelian local field theory at the shorter scale than R M - 
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FIGURES 

FIG. 1. (a) The dual Wilson loop {W D (R,T)) M A as the function of its area R x T (b) 
(Wd(R,T))ma as the function of its perimeter L = 2{R + T) in the monopole part in the MA 
gauge in the SU(2) lattice QCD with 20 4 and (3 = 2.3. For large loops as R,T > 3, (W D (R,T)} M A 
seems to obey the perimeter law rather than the area law. 

FIG. 2. (a) The inter-monopole potential Vm{t) as the function of the inter-monopole distance r 
in the monopole part in the MA gauge in the SU(2) lattice QCD with 20 4 lattice and (3 = 2.2 ~ 2.3. 
For comparison, we plot also the linear part of the inter-quark potential V£ mear (r) = or with 
a = l.OGeV/fm by the straight line, (b) The detail of the lattice QCD data for the inter-monopole 
potential Vm{t). 

FIG. 3. The analysis for the shape of the inter-monopole potential Vm{t) in the SU(2) lattice 
QCD with 20 4 and (3 = 2.2 ~ 2.3. The solid curve denotes the simple Yukawa potential Vy(r) with 
the dual gluon mass = 0.5 GeV. The dotted curve denotes the Yukawa-type potential V(r; Rm) 
including the magnetic-size effect. The lattice data of the inter-monopole potential Vm{t) seem to 
be fitted by V(r; Rm) with the effective monopole size Rm — 0.2fm in the whole region of r. 
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